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Abstract 

Results for four-, five-, and six-parton tree amplitudes for massive quark-antiquark scattering 
with gluons are calculated using the recursion relations of Britto, Cachazo, Feng, and Witten. The 
required diagrams are generated using shifts of the momenta of a pair of massless legs to complex 
values. Checks verifying the calculations are described, and a simple formula for the shifted spinors 
of an internal gluon is presented. 
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I. INTRODUCTION 



Stimulated by Witten's introduction of a string theory on twistor space dual to pertur- 
bative gauge theory [l], significant recent progress has been made in the calculation and 
understanding of the structure of gauge-theory amplitudes. By revealing the maximally 
helicity violating amplitudes to be primitive vertices from which an amplitude can be com- 
puted, Feynman diagram calculations can be replaced by the MHV rules of Cachazo, Svrcek, 
and Witten relating the desired tree amplitude with products of propagators and on-shell 
MHV amplitudes that possess fewer legs. Britto, Cachazo, Feng, and Witten (BCFW) s], 4| 
later proved that a tree amplitude is given by the sum of products of a propagator and two 
on-shell amplitudes with fewer legs and shifted, complex momenta. This recursion yields 
very compact expressions for an amplitude. Recursive methods for the calculation of leading- 
order QCD processes have existed for many years 0, B], initially formulated by Berends and 
Giele. In contrast to the Berends-Giele recursion relation which is based on off-shell vertices, 

the building blocks for the new methods are on-shell tree amplitudes. On-shell methods have 

Q 

been previously used at loop level via the unitarity method [7|. 

The BCFW recursion, while initially derived for pure gauge boson tree amplitudes, has 
been generalized to include massive particles with spin through the work of Badger, Glover, 



Khoze, and Svrcek [8[ . Thus the new recursive methods can be brought to bear on a variety 
of phenomenologically interesting problems to be encountered at the Large Hadron Collider 
such as top quark, vector boson, and possible supersymmetric processes. The multigluon 
scattering with a single external massive gauge boson or Higgs boson has been calculated 
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121] as well as processes with a pair of massive scalar particles 
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supersymmetry Ward identities of supersymmetric QCD have been 
interactions, yielding th e q uark am plit udes for certain helicities 



gauge theory 
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35, 



36 



23, 
28, 



24 
29 



2 

301 



31 



_applied to these scalar 
14| . Besides tree-level 



26[ | , th e BCFW recursion 



32l . l33l . [SJ] and gravity 



The present paper is concerned with the tree-level scattering of a massive fermion pair 
with massless gauge bosons. Results for qq ggg were presented by Ozeren and Stirling, 
but these authors report that difficulties are encountered with the BCFW recursion so that 
Feynman diagrams were necessary when gluons are exchanged for certain helicity and spin 
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configurations 2l|. Here we remedy this and present an explicit formula for the spinors 
associated with the massless gluon exchanged within a diagram. In this paper, the massive 
quark-gluon tree amplitudes with six and fewer partons are presented. The necessary re- 



cursion diagrams are found solely via the BCFW method with spinor shi 



'ts on two of the 



external massless gluon legs. Recent work by Schwinn and Weinzierl 40| also deals with 
scattering of massive quarks which also avoid the potential difficulties. In that reference a 
formula is given for the scattering of a pair of massive quarks with an arbitrary number of 
gluons in a specific helicity configuration, one negative-helicity gluon and the rest positive. 
Here we provide all the remaining helicity configurations up to amplitudes with four external 
gluons. 

II. REVIEW OF THE SPINOR-HELICITY FORMALISM 

The tree amplitudes presented here are the color-ordered partial amplitudes An, con- 
taining the kinematic data, from which the full amplitude M„ with color information is 
determined by the color decomposition 4^ . 

The notation A;^* labels the k-th particle with spin A^, Sn-2 is the group of permutations 
on n — 2 labels, and the quark colors are ^2 and ji for the quark labelled 2 and antiquark 1. 
The color generators are normalized as tr T°-T^ = 5°"^. 

In the spinor helicity formalism, spinors for the particle with momentum pi are denoted 

u'{Pi) = u-{pi) = {i\ . (2.2) 

For massive spinors, the same angle and square brackets denote spin with respect to a fixed 
axis rather than helicity. These massive spinors have four components and satisfy the Dirac 
equation, and they are interpreted as Weyl spinors only in the massless limit. 
Spinor products take the form 

u-{PiW{Pj) = {ij), u+{pi)u-{pj) = [ij], 

u+{pi)u+{pj) = [ij), u~{pi)u^{pj) = {ij]. (2.3) 



We use the formalism of Kleiss and Stirling [41|] for the construction of massive spinors, 
which yield 



iPj ■ ko){Pi ■ h) - {Pi ■ ko){Pj ■ h) - if^fiupaKPiPj^i 

jPi ■ ko){pj ■ ki) - {pj ■ ko){pi ■ ki) - ie^^p^k^p'^p'^jk^ 
ViPi ■ ^o){Pj ■ ko) 



^■js^ ^ ^ijPj ■ ko)+mj{pi ■ kp) ^^^^ 

V(Pi ■ ^o)(Pj ■ ^o) 

where ko and ki are vectors satisfying k^ = 0, kf = —1, and k^ ■ ki = 0. The vector ko 
corresponds to the axis of spin quantization for the massive fermions. Imphcit in these 
spinor-product formulas is their analytic continuation under crossing symmetry. A past- 
pointing momentum pi is replaced with —pi while each formula acquires an overall factor of 
i, and rrii becomes — mj (m > for all legs). 

The final spinor product required is whose formula is determined by the spin-sum 

rule for a particle with momentum p, \p)\p\ + \p] {p\ = ^ + rrip, so that 

(z|^|j] = {ip)[pj] + {ip]{pj] - mp{ij]. (2.5) 

If p corresponds to an antiparticle, nip i— —rrip. The shorthand is used for this spinor 

product. Under complex conjugation, the spinor products transform as 

(ab)* = [ba], (2.6) 
[ab)* = {ab], 
{a\k\b]* = [a\k\b). 

The Schouten identity for two-component spinors is 

\a) {be) + |c) {ab) + \b) {ca) = 0, (2.7) 

and we have spinor identities for massive legs i and j, 

[ab]({ia){bj] - {tb){aj]) + {ab) ( {ta][bj] - {zb][aj]) = 2pa ■ Pb{tj], (2.8) 



mi{bc) + [ic){ib) + [ib){ci) = 0, (2.9) 

and 



[ia)[jb) = [ib)[ja). 
4 



(2.10) 



Complex conjugation of these identities yields their square-bracket counterparts. 

For helicity-labelling conventions we consider all the gluons to be outgoing so that a 
helicity label on a gluon denotes its helicity leaving a vertex. The massive quarks are both 
considered to be incoming. Thus momentum conservation means pi + P2 = Ps + ■ ■ ■ + Pn- 
Our fermion spin labels are therefore opposite to the helicity labels for outgoing particles in 
the massless limit. 

Simple relations exist between n-point amplitudes with different helicities and spins. 
Complex conjugation of an amplitude flips the spin of every particle, 

(igS 2g^, 3^^, . . . , n^") = y4„ (ig 2g ^^^5 3g . . . , rig , (2.11) 

and the spin of an individual quark may be flipped by replacing bracket spinors with angle 
spinors, or vice versa, so we have 

^n(i^2^^3^^...,n^'^)|[l|^^l| = A„(l;^,2^^3^^...,n^"), 

A„(l^^2^,3^^...,r^^")||2]^|2) = A„(l^^2;^,3^^...,n^"). (2.12) 

|2) ^ |2] 

Thus, for example, complex conjugation of the amplitude An (ig , 2^ , 3*^, . . . , n^") followed 
by the replacements (1| ^ [1\ and |2) i-^ |2] flips the helicity of only the gluons, yielding the 
amplitude An (l^, 2+, 3""*^, . . . ,n~'^"). Explicitly following the terms [la] and [la) through 
this composition of maps we have 

[la] ^ -(la) -[la), (2.13) 

[la) A (la] ^'J^'' [la]. 

Every possible combination of spins can be obtained from the amplitudes 
An (iq ,2^ ,3g^ , . . . ,ng") by a combination of complex conjugation and quark spinor 
replacements. Thus in presenting the results for the amplitudes we need only present 
the lg,2g amplitudes and half of the possible gluon helicity configurations. The spinor 
identities (12.71) . (12. 8p . (12.91) . and (I2.10p can only be used when the massive fermions have 
particular spins, so their use will invalidate the above relations between amplitudes. We 
present the amplitudes in a form where no spin-dependent identities are applied in order to 
allow relations (12. lip and (12.12P to be used. 
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The amplitudes presented in ref. 40| use a different spinor formalism which must be 
accounted for in comparison of the amplitudes. The spin quantization axis ko is chosen to 
be one of the massless gluon legs, labelled as leg q, and massive vectors are projected to the 
light cone with this leg q, 

Massive spinors are defined as 

«-(*) = 1*1+ (^k). (2-15) 
which leads to the spinor products 

{»} = fol = 

III. THE BCFW RECURSION FORMULA 



The BCFW tree-level recursion formula follows from the basic complex analytic proper- 
ties of amplitudes in gauge theory. Consider a tree-level color-ordered gauge-theory ampli- 
tude An{pi, . . . ,Pn) with complex momenta obtained by shifting a pair of spinors from the 
massless external legs i and j, 

^=\^) + z\J), \j] = \j]-z\t]. (3.1) 

Under this spinor shift, the momenta shift as 

Pi = K)N 

P. = |J)[J|-^|J)[^I- (3.2) 

To be concise we denote this choice of shifted legs "(ij]", but from context this should not 
be confused with the spinor product of eq. ( 12. 3p . 

This BCFW shift keeps all momenta on shell, and total momentum is still conserved. 

Pi +Pj = Pi +Pj- 
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FIG. 1: An illustration of the recursion formula for cases where the shifted external legs pi and pj 
are chosen to be adjacent with i < j. The sum over the internal particle's helicity is implicit. We 
denote the BCFW diagrams, from left to right, as Vi, 'D2, ... 

The amplitude is now a meromorphic function An{z) with simple poles in z when the sum 
of adjacent particles' momenta, denoted P, is z-dependent and on shell, 



= 



m 



- + z{j\P\i] 



(3.3) 



Provided that An{z) vanishes as z —* 00, we have J A{z)/z = 0. Applying Cauchy's theorem 
to J A{z)/z yields residues for z = (the desired unshifted amplitude) and the values of z 
w? = 0. The factorization of amplitudes at these poles yields the BCFW 



for which P^ 



formula 



1 



P^ — nip 



X Ar\ 



right 



(-p-^...,p«), (3.4) 



partitions h 

as illustrated in fig. [H Each term appearing in the recursion formula is one of the residues of 
A{z) I z. The momentum P is the shifted total momentum leaving from the right amplitude's 
external legs, and the recursion formula sums over the helicity of the internal state with 
momentum P. The sub-amplitudes v4ieft/right are to be evaluated with the shifted spinors at 
the value of z for which P^ — rr? = 0. 

Valid recursion relations are obtained only for certain helicities of the shifted external 



legs, where An{z) vanishes as z —>■ 00. Ref. [40j derives the valid choices of legs to shift 
for both massive and massless particles. We shift only massless gluons in this paper, so it 
suffices to note that the recursion is valid for the "(ij]" shift in eq. (13. ip provided that the 
shifted legs i and j have helicities 



{hi, hj 



:+,+), (+,-: 



or 



(3.5) 
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By applying the spin-sum rule to combine the spin states of an internal fermion, rather 
than considering spin states separately, a term with a fermionic-particle pole in the recursion 
formula can instead be written as 

where P* denotes that the internal spinors have been stripped off the amplitudes. We use 
this method, from ref. L 



to calculate the BCFW diagrams where an internal fermion is 
exchanged between trees. An amplitude which is to be fed into the recursion should not 
be simplified in advance using any identities which depend on the spin of the massive legs. 
Doing so violates the rule that both spin states of an internal fermion should be summed 
over, not just the spin for which a particular identity holds. 

Consider the spinors | P) and | P] which appear in the recursion diagrams for the exchange 
of a massless particle. Because the internal particle appears with opposite helicity at each 
sub-amplitude, the overall amplitude will have degree zero in the P spinors. Thus identities 
such as {aP)[Pb] = {a\P\b] may be used to eliminate the individual spinors. However, to 
facilitate computer programmable calculations of the diagrams, it is convenient to have a 
direct substitution for \P) and \P]. Moreover, an explicit formula for the P spinors solves 



the problem of dealing with massive spinor products such as [IP) encountered in ref. |21 |. 
To achieve this, we use the two-dimensionality of massless spinors to write the P spinors in 
terms of the shifted legs i and j (the shift is "(zj]"), 

\P) = a\t)+P\j), \P]=^\t] + S\j]. (3.7) 

Then by using the identities 

{P\P\t] = {j\P\P]=0, 

{j\P\t] = {j\P\t], (3.8) 

we find 

^ iP\^])i{J\P), (3.9) 



{j\P\^ 

providing a simple formula for the shifted spinors of an internal gluon. The overall prefactor 
may be associated with either \P) or \P]. This expression for P is valid no matter on which 
side, left or right, the legs i and j lie. 
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FIG. 2: The single BCFW diagram contributing to all the four-point amplitudes. 
IV. qqgg AMPLITUDES 

As input for the single BCFW diagram of fig. [2] contributing in the amplitudes with 
two gluons we use the polarization vectors of positive- and negative-helicity gluons with 
momentum p, 

fiP.Q) = + (4-1) 

with the three- vertex factor The spinors q refer to an arbitrary null vector. Different 
choices of reference spinors lead to gauge-equivalent polarization vectors. 
The four-point amplitudes with opposite-helicity gluons are 

^,(I-.2-.3j.4;)H4|l|3lMM±MM (4.2) 

and 

^,(I-,2-,3,^4j)^|4|l|3)»MM. (4.3) 

^ 4p3 • p4Pi ■ P4 

The four-point amplitudes with identical-helicity gluons are given by 

.4I,^2,^3,,4,)^™(43)M|™ (4.4) 

and 

.41^,2,3^4;)., ,34| l'^>!g™ . (4.5, 

Using eq. (12.101) this amplitude vanishes but only for this specific choice of quark spins. To 
obtain the amplitudes with the other spins, one starts from eq. fl4.5p and applies eqs. fl2.1ip 
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FIG. 3: The BCFW recursion diagrams required for the five-point amphtudes under the shifts 
"(45]" or "(54]". These diagrams are relevant for the gluon hehcities ( ) and (H ). 

and fl2.12p . One must include this amplitude in this form as input to higher-point recursion 
diagrams with internal quarks, since the internal quark is summed over all spins and not 
only the spin for which this amplitude happens to vanish. 

Likewise note that the Schouten identity could be apphed in eq. (14. 4p to yield [13] [42] — 
[14] [32] = [12] [43], thus cancelling the unphysical double pole [34]^ in the denominator. 
However this identity can be applied only for the particular quark spins chosen, so we leave 
the expression as is to avoid invalidating the next level of recursion. 



V. qqggg AMPLITUDES 



We now write the five-point amplitudes beginning with As (l^ , 2^ , 3^ , 4^ , 5^ ) calculated 
from the shift choice "(45]", with the result 



^5 (iq ' 2g , 3g , 4^ , 5g ) 



'"((34>-.|fjf(35> 

m[15)[43][42] - [14][3|1|5)[42] + [14][4|1|5)[32]' 
2p2-P3[34]2[45]2pi-p5 



(5.1) 



With the shift "(45]" we find 



^5 (l-,2-,3+,4^-,5,-) 



[3|2|4) + ^^,Vif' [3|2|5) 



X 



+ 



[4|1|5) 

(-[14][3|l|5)+m[15)[43])((42] + te.(52 



(5.2) 



(- [15)2p2-P3 + m[14](54))[32]" 



8P2 -^3(^3 -^4 + fi||[4|3|5))pi •P5[45] 



m(45)' 



(43)(pi+p2)^(2pi-p5 + ||(5|l|4]) 



[1P^][52]-[15][P«2] 
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FIG. 4: The BCFW recursion diagrams required for the five-point amphtudes under the shifts 
"(34]" or "(43]". These diagrams are relevant for the gluon hehcities ( — I — ) and ( h). 



For the second term, which arises from the first BCFW diagram in fig. [3l we have z = — -jHy 



so that 



' ^ ' ^ (35)' ' (35)' 
|5]H5] + |||4], 

Pr = P3+P4=|3)[3| + |4)[4| = |3)([3| + ||[4| 



(5.3) 



where we have used the Schouten identity. Thus, for example, [oPr] = [a3] + |||y [a4] . We use 
this notation throughout the paper, indicating the spinors to be substituted for the internal 
massless momenta P or Pr by explicitly writing /' as a bi-spinor. 
The shift "(43]" leads to 



A, (1,^,2,-, 3;, 4+ 5;) = [4|1|5) 



m[14](53)[42] + [15) [42] [4|2|3) - [14](32][4|1|5) 



+ 



m ■ P5P2 ■ P3 [P4 ■ P5 I 

m(35)^ 



£^[4|5|3) [34] 



(5.4) 



(45)(34)(pi +p2)^(2p2 ■ P3 + ||(3|2|4]) 



[13][P2]-[1P][32] 



where, in the second term, which arises from the first BCFW diagram in fig. |H z 
which means 



■ (35) ^ 



|3]H3] + |||4], 



P = P4+P5 = |5)(^|^[4| + [5| 



(34) 



(5.5) 
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The amplitude (l^ , 2^ , 3^ , 4^ , 5+) is calculated from the shift "(43]' 
A (i-,2-,3;,4;,5+) = ((4|1|5] + ^^(3|1|5] 



(3|2|4] 

[15](m(43)[42] + 2pi -^5(32]) + ([14) + gm[13)) ((3|2|5] [42] + m[54](32]) 
8pi ■ P,P2 ■ P3(P4 ■ P5 + fgj(3|5|4])[34] 
m(34)3 



(5.6) 



(45)(pi + p2r(2p2-P3 + ||(3|2|4]) 



X 



[13][P2]-[1P][32] 



where the shifted spinors in the second term are the same as in eq. (15.51) . 



VI. qqgggg AMPLITUDES 



The six-point amplitudes with zero- and one-positive-helicity gluon are presented first. 
We have 



^6 (V. 2,-, 3,-, 4;, 5;, 6;) = -m((34) + ^^(35) 



1 



X 



([16)2P-P5 +m[15](65))m[43][42] 



2p2-P3[43]2[45]2P-p5[65]2pi-p6 

(6.1) 



+ (-[15][4|1|6) + m[16)[54])([3|P|5)[42] - [4|P|5)[32]) 



where the shift "(56]" is used, so that 



15) - 15) + ^16), 16] H6]- 1^15], 



[5|1|6) 

[a|P|5) = [a|2-3-4|5), 



[5|1|6) 



-P2 ■P3-P2-Pi+P3-P4- 



(6.2) 



We also use the "(56]" shift to obtain Aq (l^ , 2^ , 3+, 4^ , 5^ , 6^ ) . This gives multiple non- 
vanishing diagrams as shown in fig. [5] which we denote as 



where 

vt— = 



-m- 



A, (1-, 2-, 3+ 4;, 5;, 6;) = P+— + Vt— + Vt"-, 

■- 0, and the second diagram under the "(56]" shift yields 
(45)' 



(6.3) 



(5|P|6](3|P|6](34)p|4_52pi-pg 



[1|P^|6)[62] + [16](6|P^|2] 



(6.4) 
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FIG. 5: The BCFW recursion diagrams required for the six-point amphtudes under the shifts 

"(56]" or "(65]". This set of diagrams is relevant for the amphtudes with gluon hehcities ( ), 

(+ ),(- + --),(- + - + ), and (- + + - ). 



where z 



Pa, 4, 5 



(6|3+4|5] 



SO that 



|5) 
(a|P|6] 



|5)- 



^3,4,5 



■|6), |6] = |6] 



^3,4,5 



(6|3 + 4|5] 
H(l + 2)|6], 

PR\h\ = (63)[36] + (64)[46] + (65)[56]. 



(6|3 + 4|5] 



|5], 



For the third diagram, , we find 

vt — - 



[3|2|4) + ^^^[3|2|5) 



[4|P|5) 



1 



X ■ 



2P2 ■ P3 ( 2p3 ■ PA + ^[4|3|5) ) 2P ■ P5 [54]2pi ■ p6[65] 



(6.5) 



(6.6) 



[3|P|5) (-[15][4|1|6) +m[16)[54]) +m[43]([16)2P-p5 +m[15](65))) 

xf(42] + ^^^(52]) 
V [4|P|5) ^ 

- ([16)2P-pg +m[15](65))2p2-P3[32] +m(-[4|l|6)[15] +m[54][16)) (54) [32] 

(45)3 



+ m 



[65]2pi ■p6(43)(pi +P2 -P6)'(2Pi -Ps - 2P5 -Pe + ||((5|1|4] - (5|6|4])) 



X 



(m[16)([53] + j||[54]) - [15]((6|1|3] + ^(6|1|4])) ([52] + ||[42]) 
(35) (35) (35) 

- (m[16)M[54] - [15]((61115] + M(6lll4])) ([32] + i§|[42]) 
(35) (35) (35) 
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In Vt we have, 



p . = -p^.p^-p^.p^ + p^. p^^ 
[a\P\5) = [a|2-3-4|5). 

From the shift "(56]" we have 



^6 (l-,2^-,3,-,4+,5,-,6,-)=Pr+- 
where the first diagram in fig. O is 

{56>A\, (34) + 1^(36) 



—m 



(46) 



X 



(43> J (64> 

m[16)[PR3][P^2] - [1P^][3|1|6)[P^2] + [1P^][P^,|1|6)[32] 



with |Pr] = |4] + -ply 1 5]. The second diagram yields 



(35)^ 



(34)(45)(5|P|6](3|P|6]pi,4,52pi-pg 
and the third diagram 

[4|P|5) 



X 



[1|P^|6)[62] + [16](6|P«|2] 



+ 



2P ■ pg2p2 ■ P3 (2p4 ■ PS + |ili[4|5|3)) [34]2pi ■ p,[65] 
;-[15][4|l|6) +m[16)[54]) ([4|P|5)(32] + m(53)[42]) 



+ ([16)2P-pg +m[15](65))[4|2|3)[42] 

m(35)3 

(54)(34)(pi -P6)'(2P2 -Ps + |g(3|2|4])2pi ■p6[65] 

( - [15]((6|1|3] + ^(6|1|4]) + m[16)([53] + ^[54])) ([52] + M[42]) 
(35) (35) (35) 

- ( - [15]((6|1|5] + ||(6|1|4]) + m[16)||[54]) ([32] + ^[42]) 
(35) (35) (35) 



where the spinor shifts in V2 ^ and "Dg 
tively. 



are the same as eqs. fl6.5p and (16. 7p . 
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Choosing the shift "(43]" leads to the diagrams in fig. [HI 

A, (1-, 2-, 3;, 4;, 5+ 6;) = V"+- + ©2^-+- + VI 

where we find 

[5|1|6) 



— +- 



2pi ■ p,2'P ■ P2 [2p, ■ p, + |^[5|6r4) j [45]2p2 ■ PsM 
([16)[5|P|4) + m[15](64)) ([5|2|3)[42] + m[54](32]) 



X 



+ m 



- [15][5|1|6)(2P-P2(32] +m(43)[42]) 

(46)4 



X 



(65)(45)(pi + j)2 -P3)'(2^ -Pi + |iy(4|i'|5])[34]2p2 -Pa 
[14] + ||[15]) (((3|2|6] + ||(3|2|5])[42] + m([64] + ||[54])(32] 

- ([^6] + Ht^^O ((^^1^1^] + g(3|2|5])[42] + ^||(32] 



(46) 

The shifts for are 



(46) 



|4) = |4) + ^|3), |3] = |3]-^|4], 



[4|2|3) 
[5|P|4) = [5|6-1|4), 



[4|2|3) 



The second diagram, V2 ^ , is 

(46)4 



V. 



— +- 



— m- 



(45)(56)(6|P|3](4|P|3]2p2-P3P4, 



5,6 



[13](3|P|2] + [1|P|3)[32] 



where D2 has z 



Pi.5,6 



(3|5+6|4] 



and 



(3|P|6] = (34) [46] + (35) [56] + (36) [66]. 



For the third diagram we have 

2pi-P6 



1^3 = —m 



2P2-PsP3,m[-P612pi-P6{54> (35) 



^3,4,5 I 



m[16)^[P2] - [1P](6|1|31[P2] + [1P](6|1|P] [32] 
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FIG. 6: The BCFW recursion diagrams used for the six-point amplitudes under the shifts "(34]" 

or "(43]". These diagrams are apphed for the amphtudes with gluon hehcities ( 1- — ), ( \-), 

and (-- + + ). 

where in the third diagram z = — -jHy and \P] = |5] + |||y|4]. 

The "(43]" shift yields the final single-positive-helicity gluon amplitude at six points, 

A (i-, 2-, 3;, 4;, 5;, 6+) = V^-+ + V,-^ + P3— +, (6.18) 

where ^ vanishes. For the remaining diagrams, V'l ^ and V2 we have 



1^1 



+ 



[6|1|5) + 



^^[61114) 

[5|P|4> L^l-"!^/ 



X 



2pi-p62P-P4 (2P5-P6 + |^[5|6|4)) [45][34] (-2^2-^3) 
[16]2pi ■ p6 + m([15) + ^^[14)) [65]) (2P ■ (32] + m(43) [42] 

[15) + [14)) [6|P|4) + m[16](54)) f [5|2|3)[42] + m[54](32] 



(6.19) 



[5|P|4) 



+ m 



(45)^ 



(65)(pi +P2 -P3)'(2^ - Pi + |iy(4|i^|5])[34]2p2 - Pa 



X 



(46) 



(46) 



(46) 



and the second diagram 



—m- 



(45)= 



[13](3|P|2] + [1|P|3)[32] 



(6.20) 



(56)(6|P|3](4|P|3]2p2-P3Pl5,6 
The shifted spinors in the first and second diagrams, ( 16.19^ and ( 16.20p . are identical to those 
found in Aq (1~, 2~, 3~, 4~, 5^ , 6~), eqs. f l6.14p and fl6.16p respectively. 
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Finally we write the amplitudes with two negative-helicity gluons, beginning with the 
amplitude 



^6(1-, 2- , 3; , 4; , 5+ , 6+) = vr^+ + + n 

where the "(34]" shift is chosen. The first diagram, "Df is 

(5|1|6] 



++ 

3 ' 



-++ 



(4|P|5] + |^(4|P|6] 



X 



2P • P4 (2p3 • p5 + |frg(5|4|6]) 2pi • pe{65) (-2p2 • P3) [34] 
[15)(4|1|6] +m[16](54)' 



X ( (3|2|5][42] + m[54](32] + ((3|2|6][42] + m[64](32]) 



+ 



(-[16]2P-pj +m[15)[65]) (^2P • (32] + m(43)[42]) 

m[56]'^ 

(pi+n-Ps)^ {iPi-Pe + |}<5|l|6l) [54]143]2p2 .p3 
For the diagram Vi~~^~^ the spinors are shifted with 



X 



z — 



-2P2 • P3 

(4|2|3] ' 



(4|P|5] = (4|6-1|5], (a|P|4] = (a|5 + 6- 1|4], 

P-PZ = Pi • P5 + Pi • P6 - P5 • P6- 



(6.21) 



(6.22) 



(6.23) 



We note that the second term in T>i arises from a vanishing term in 

(l~, 2~, 3~, 4+, 5+). This five-point term is zero for the special quark spins chosen but 
contributes a residue at six points since the recursion sums over both spin states of an 
internal quark. The second diagram yields 

(3|2|5](54) + (3|2|6](64) [65] 



X 



(pi+P2)'2p2-P3Pl,5,6 (45) (56) [54] 
[13) ((45)[52] + (46)[62]) - ([15](54) + [16](64)) (32] 



(6.24) 
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where the P spinors have been ehminated using {aP)[Pb] = (a|4 + 5 + 6|6], and z = j^^!^^' 
For the third diagram we have 



T)—++ 



(3|2|5] + |fg(3|2|6] 



[35] 



2P2 -Pa {^P-P3+ f[S(^|3|6]) 2pi -^6(6^) [45][34] 



(6.25) 



X 



-[lP)(3|l|6]+m[16](P3)) ([52] + ^|^[62]) - {[16]2p, ■ p^, - m[lP)[65]) (32] 



where z 



45 



35 



(P|l|6]' 

By using the Schouten identity we find 



[34] 
[35] 



|4) + |5) [5 



(6.26) 



The next amphtude is 



A, (i-, 2^-, 3^-, 4+, 5;, 6+) = Pr+-+ + ^2-+-+ + V, 



(6.27) 



where the shift "(65]" is chosen to find 

(3|2|4] + |Sl^3|2|6] 



■+-+ 



[A6f 



X 



2P2 • PS {2Pn ■ PS + ^^{PRim) 2pi • PePh,e [^^^t^^^ 

(-[1P«)(3|1|6] + m[16](P^3)) ([42] + ^^[62]) 
^ ^ ^ {Pr\1\Q\ ^ 

(-[16]2p2-P3 + m[lP,^)[64]) (32] 



(6.28) 



+ 



where z = ||||, so then Pr = (|4) + 1^15) j [4| and 2Pr-p^ ^ ^'3 4 5- -^^^ second diagram, 
V2~^~~^, we have 



(35)^ 



X 



(34)(45)(5|3 + 4|6](3|4 + 5|6] p^,^, (p, + p,f 2p, ■ 



(6.29) 



x[6P](P|2|6] X 



[1P)[P6][62] - [16][6P](P2] 
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where z = (^^\^'^'1\q-\ and |P)[P6] = — {ps + P4: + P5) |6]. The third diagram is 



T,-+- + 
^3 



(5|6-1|4] 



2P ■ P5 2p2 ■ P3 (2p4 ■ Ps + Tig(3|5|4]) [34] {-2p, ■ p,) (56) 



(6.30) 



— m 



([16](5|P|4] +m[15)[64]) (m{53)[42] + (5|P|4](32]) + [15)(5|P|6](3|2|4][42] 

W 

(56)2pi ■p6(54)(34)(pi -P6)'(2P2 -Ps + S(3|2|4]) 



X 



(m[15)([63] + ||[64]) + [16]((5|6-l|3] 



(35> 

||(5|6-1|4])([52] + ||[42]) 



- ("^[15)([65] + ||[64]) + [16](2P -pg + ||(5|6 - 1|4])) ([32] + ||[42]) 

where z = -^gyf • Then P = p^-pi so that (5[P[6] = (5[6- 1[6] and 2P ■ p^ = -2p2-p3- 
2p2 ■ Pa + 2p-i ■ Pi. 

The final six-point amphtude with two negative-hehcity gluons is found from the shift 
choice "(56]", 



A, (1-, 2-, 3;, 4+, 5+, 6;) = V^++- + I?2 



(6.31) 



We find for the first diagram in fig. [5] 



(46) 



X 



(6[1|P, 



R\ 



X 



(56)(45) 2p,.p^,2p,.p, (2P« -pg + ||g|(3[6[P^]) [3P«] 



m[lPij](63)[PR2] + [16)(3[2[P«][P^2] - [lPR](6[l[Pi^](32] 

(36)^ 



(6.32) 



(34)(45)(56)(pi + p2r(2p2 -^3 + ||(3|2[Pr]) 



X 



(1131 -H Iflli^HlX^S] + |||P„21) - (116] + |||iPh1)(|321 + M 



I^'r21) 



where ^ = — -jlgy and Pr = [4) ( [4[ + ||gy[5[ ). The second diagram is 



n 



-++- 



[45]= 



X 



(6|P|5](6|P|3][34] p|4^52P-pg2pi-pg 



(6.33) 



x(6|2|P](P6) X 



[1P](P6)(62] + [16)(6P)[P2] 
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where z = - (6^3+^5] > -'^P ■ H = (Pi + P^f , and |P](P6) = - (ps +P4 + P5) |6). The third 
diagram, Dg"*"^", is 

(3|2|4] + ^(3|2|5] 

= ^ ^^^^4 (6.34) 



2p2 ■ P3 ( 2p3 ■ P4 + ^(4|3|5] ) 2P ■ (54) [65]2pi ■ p. 



X 



[16)(4|P|5] I (3|P|5]( [42] + ^^^[52] ) +m[54](32] 



X 



(4|P|5] 

+ [15] (^(6|P|5] (^m(43) ([42] + ^^[52]) - 2p, ■ ps{32] 

+ m(64) (^(3|P|5](^[42] + ^^[52]^ +m[54](32] 

m[45]^ 

{Pi +P2- n)' (-2^ ■ P5 - {i(4|^|5]) [43] [56]2pi ■ p, 
[16) (^((3|P|5] + ||(4|P|5]) ((52] + M(42] 

- (2P . + ||(4|P|5]) ((32] + ^(42])) + -[15]||(62] 
where z = ^f^, 2P ■ p^ = {p2 - ps - Pif , and (a|P|5] = (a|6 - 1|5]. 

VII. CONCLUSIONS 



We used the BCFW formahsm to calculate the scattering amplitudes with a pair of 
massive fermions and four or fewer gluons by shifting adjacent pairs of massless gluons. 
All the helicity configurations are computed in this way. By deriving an explicit formula 



n2 



2l| with the 



for the shifted spinors of an internal gluon, no subtleties are encountered 
massive-spinor inner products. 

We have performed a number of checks on our results. By performing a variety of shifts 
we numerically confirmed that the residues corresponding to all possible channels are correct. 
We note that shifts where the amplitude is not well behaved for large z can be used for the 
purpose of checking a putative amplitude, even though an on-shell recursion constructed 
from such a shift would drop terms coming from the additional large z contribution. To 
perform the check we shift a different pair of spinors than the ones used to construct the 
amplitude via on-shell recursion and extract each residue which we then compare to the 
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corresponding BCFW diagram. This verifies that the answers are consistent, so that we 
obtain the proper residues no matter which pair of external gluons are shifted. This check 
is equivalent to checking all factorization limits of the amplitude. 



We have also checked that our results agree with the amplitudes presented in ref. 
where formulas for the amplitudes with identical or all but one identical helicity gluons were 
derived from massive scalar amplitudes via supersymmetric Ward identities. The five-point 



amplitudes match the results of ref. j2lj, which were calculated using BCFW recursion and 
Feynman diagrams for certain amplitudes. 

The amplitudes presented in this paper are a subset of those needed for heavy quark 
plus two-jet production at the LHC The calculations of this paper confirm that the BCFW 
recursion relations are an effective way to compute amplitudes with massive quarks. 
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